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ABSTRACT 

Imaging data from the Sloan Digital Sky Survey is used to measure the empirical size-richness 
relation for a large sample of galaxy clusters. Using population subtraction methods, we determine 
the radius at which the cluster galaxy number density is 200O~^ times the mean galaxy density, 
without assuming a model for the radial distribution of galaxies in clusters. If these galaxies are 
unbiased on Mpc scales, this galaxy-density-based i?200 reflects the characteristic radii of clusters. We 
measure the scaling of this characteristic radius with richness over an order of magnitude in cluster 
richness, from rich clusters to poor groups. We use this information to examine the radial profiles 
of galaxies in clusters as a function of cluster richness, finding that the concentration of the galaxy 
distribution decreases with richness and is systematically lower than the concentrations measured for 
dark matter profiles in N-body simulations. Using these scaled radii, we investigate the behavior of 
the cluster luminosity function, and find that it is well matched by a Schechter function for galaxies 
brighter than Mr = — 18 only after the central galaxy has been removed. We find that the luminosity 
function varies with richness and with distance from the cluster center, underscoring the importance 
of using an aperture that scales with cluster mass to compare physically equivalent regions of these 
different systems. We note that the lowest richness systems in our catalog have properties consistent 
with those expected of the earliest-forming halos; our cluster-finding algorithm, in addition to reliably 
finding clusters, may be efficient at finding fossil groups. 
Subject headings: galaxies: clusters: general — cosmology: observations 



1. INTRODUCTION 

Galaxy groups and clusters appear in many guises. Ob- 
servationally, they can be identified as pools of X-ray 
emitting gas, collections of galaxies, Sunyaev-Zel'dovich 
decrements of the CMB, or strong features in the gravi- 
tational shear field. Theoretically, they are identified as 
the largest over-dense 'halos' of dark matter. Mass is the 
key defining attribute of a galaxy cluster. The evolution 
of the cluster mass function and its variance plays an im- 
portant role in constraining cosmological pa rameters de- 
scrib ing dark energy, such as w and JIa fe.g-. lLima fc Hu.l 
[2OO5I), and large-scale structure, such as Vim and ag (e.g. 
^ahcall ct al. 2003). Mass estimates and studies of clus- 
ter members both rely upon knowing the size of a cluster 
of given mass. To define the mass of a cluster, typically 
a cluster radius is specified through some prescription, 
and the total mass taken to be the mass contained within 
that radius. Various operational definitions are used to 
determine cluster size and mass, but all of these show 
a regular increase in cluster size with mass. In order to 
study objects in a wide range of masses in a meaningful 
way, we must be able to determine an appropriate size 
scale for these clusters and groups; we may then use this 
characteristic scale as an aperture within which to make 
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comparisons. 

In numerical simulations, the precision of cluster mass 
and size measurements is limited only by resolution. 
Still, there are a variety of definitions for bo th size 
and in ass in use, as discussed in some detail bv lWhit3 
l)2001|) . One class of estimates is based on top-hat 
filtered spherical over-densities. In this model, clus- 
ters are expected to be virialized within regions where 
the enclosed mean mass densi ty exceeds the critical 
densi ty by a factor A - 200 l)Peeblesl 119931: iPeacockl 
119991) . The radius at which this over-density is reached, 
^A=200j is used as the characteristic radius of the clus- 
ter. The total mass within this radius, il/A=200j is 
used as the characteristic mass. A number of choices 
of A are in use in the literature, from an overdensity 
of 180 times the mean background (e.g. I.Tenkins et alJ 
2001; Kravtsov ct al. 2004J, 200 times the critical den- 
sity, i.e. 20007^ times the niean background (e.g. 
'Diafe rio et aI]l200H lE7 rard et al.ll2002t iKochanek et alJ 
2003 ) , to the "viri al mass" feke. Cole, FrenkI 11 Pflfit 
Bullo ck et alJl200i|) . Alternative definitions ident ify ha- 
los by 'friends-of-friends' (FOE) methods (Dav is et alJ 
Il985t Ijenkins et all 120011) . In these methods, particles 
are associated with halos to which they are linked by se- 
quences of neighboring particles. Masses for EOF halos 
are often given by the sum of member particle masses, 
but as the halos are not required to be spherical, halo 
size is less clearly defined. 

Observationally, cluster mass and size are difficult 
to measure directly; typically some mass estimator is 
adopted as a proxy, and typically a mass model is as- 
sumed to calculate the virial size. With deep observa- 
tions, Icnsing may be used to make detailed mass maps of 
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rich clusters, but this is not yet a practical technique for a 
large sample of systems spanning a wide range of masses. 
With large, shallower surveys, the lensing signals from a 
set of lens syste ms can be stacked to determine a compos- 
ite mass profile fSheldon ct al. 2001). However, the large 
number of lens systems needed for high signal to noise 
makes such a measurement difficult for examining narrow 
mass ranges. While examining the lensing signal from an 
individual cluster recovers a specific mass, that estimate 
is affected by the errors, such as those due to projection 
effects, in modeling the mass. Combining the lensing sig- 
nal from many clusters is advantageous because projec- 
tion effects are unimportant, and the stacking simplifies 
modeling of neighboring structures. While stacking does 
limit the resolution in mass for a large set of clusters, this 
limitation is due to the amount of data available rather 
than the technique used. Other observational techniques 
for estimating mass, whether using X-ray temperature, 
SZ flux, or velocity dispersion as a mass estimator, also 
rely upon either expensive spectroscopy to gain a detailed 
understanding of a few rich systems, or upon models for 
the mass distribution used to infer the total system mass 
and size. To compare observational data to theoretical 
models, and to compare the observed properties of clus- 
ters of different masses, it would be preferable to avoid 
using a model-dependent mass/radius scaling. 

With a large photometric optical survey it is now 
possible to use cluster richness to characterize galaxy 
systems without suffering from the projection effects that 
have plagued such a mass estimator in the past. Since 
the cluster catalog of, Abcll (1958), systems of galaxies 
have been sorted and compared using a variety of 
richness parameters, many of which have been based on 
the number of galaxies within a certain luminosity range 
and distance from the estimated cluster center. The 
richn ess parameters of i Abell. Corwin & Olowin ( 1 98Q.); 
ICouch et a l. (199ll: Da lton et al . (1992); Lumsdc n et al] 
( 1992[) : iLidman fc PetersonI (11996): .Postraan et al.l 
h99&): fO strander ^all (119981); lOlsen et all 119991): 
iCTladdersI (Ii200a): 'Goto et al.' (2002); 'Postman et al.' 
(1222^; iCal et alJ (2 003) arc of this type; see Bahgall 
111981 : lYee fc Lopez-Cruzi (I1999D : IBahcall et all i2(Mi 
for further discussion and comparison of some of these 
richness estimators. We can stack systems in narrow 
richness bins, and measure directly the distribution of 
galaxies in clusters over a wide range of masses. This 
galaxy distribution is used to estimate the virial size 
of these systems. Ideally, cluster members would be 
identified spectroscopically. Such data is not feasible 
to obtain for a very large sample of clusters, so we 
rely upon projected photometric data taken in multiple 
bandpasses, and correct for the foreground and back- 
ground galaxies that contaminate our line of sight to 
each cluster. 

The Sloan Digital Sky Survey (SDSS,|^rrsE^|2QO0; 
l^ouahton ct al. 2002) data offer thousands of clusters 
and groups for study, and can be used to measure excel- 
lent photometric redshifts for those objects. We use data 
from the SDSS to directly determine a size-richness re- 
lation for groups and clusters with a model-independent 
method. Since we cannot directly measure the radius at 
which the cluster has a mass over-density of AM , we in- 
stead determine the radius i?'^ at which the space density 
of cluster galaxies. A/", is over-dense by A7V. We present 



the scaling of i?200 with richness, which can be employed 
to further study galaxy clusters and their members. For 
example, to study the relationship between different mass 
estimators, such as the mass from lensing measurements 
and the total luminosity of the cluster, it is essential to 
know this scaling of cluster size with richness. 

If the distribution of galaxies in a halo traces the over- 
all dark matter distribution, our galaxy-density based 
i?2oo '^ill reflect the characteristic radius of clusters. In 
detail, the relation between the dark matter density pro- 
file and the radial distribution of a population of galaxies 
depends on a number of physical processes including dy- 
namical friction and tidal striping, and depends on the 
properties of the galaxy sample, but both simulations 
and previous observational work suggest that the distri- 
bution of galaxies in a halo at least roughly traces the 
overall dark matter distribution (e. g . , iNagai fc Kravtsovl 
l2?)0llLin. Mohr. fc Stanfordll2fl0l . 

This hypothesis is supported b y recent lensing and 
galaxy clustering measurements ()Sheldon et al.l 120041: 
Weinberg et al. 2004), which suggest that the bias of typ- 
ical SDSS galaxies is approximately one and is roughly 
scale-independent on scales larger than a few hundred 
kpc. In any case, without assuming a model for the 
radial distribution of galaxies in clusters, we can directly 
measure the radius at which the galaxy density in clus- 
ters is A times more dense than the average background. 

Our investigation does not distinguish between galaxy 
groups and clusters. Systems of galaxies come in a range 
of masses from single galaxies in larger halos up to the 
most massive of clusters. There is a clear boundary be- 
tween galaxies and systems of galaxies; there is not any 
clear dividing point between poor and rich systems of 
galaxies as observed in the optical. We use the full range 
of our cluster finder to develop our catalog of systems; 
our mass estimator spans the range from highly popu- 
lated clusters down to very sparse systems, and we refer 
to any system with two or more galaxies as a "cluster" . 
There is no strong break in the properties of the systems 
found as a function of richness. 

In [J21we describe the SDSS data used and discuss the 
cluster finding technique and richness measurement; in 
iJSlwe present and test our method of background sub- 
traction through examination of the radial profile and lu- 
minosity function. Our calculation of R^q as a function 
of richness is presented in 21 We further examine the 
radial density profile within i?2oo 33 we examine the 
luminosity function within -R200 how the luminos- 
ity function changes as a function of r/R^Q for clusters 
of different richnesses in ^ Throughout, we assume a 
flat, LCDM cosmology with Ho = lOOh km Mpc"\ 
h — 0.7, and matter density flm = 0.3. 

2. DATA 

2.1. SDSS Galaxies 

In this study, we use 39 5 deg^ of SDSS commission- 
ing data l)York et al.ll200nt) . a subset of the Early Data 
Release data. In particular, we use the 170 contiguous 
square degrees imaged September 19 and 25 1998, cov- 
ering the range 145.1° < RA < 236.0°, -1.25° < DEC 
< -fl.25° (J2000), known as SDSS stripe 10, and the 
225 contiguous square degrees imaged March 20-21 1999, 
covering the range 351° < RA < 56°, -1.25° < DEC < 
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+ 1.25° (J2000), known as stripe 82. Seeing varies on 
these two stripes from 1.0" to 2.0", and the data are 
photometrically u niform to within 3% ijHogg et alJl2001t 
iSmith et alJl2002() . Star-galaxy separation is robust to 
21.0, 21.0, 2 1.0, and 19.8 in q, r i, and z passbands 
respectively l)Scranton et al.lf2002D . which we adopt as 
the limiting apparent magnitudes for this work. All 
apparent magnitudes are measured by the photometric 
data proc e ssing pipeline using a modified version of the 
iPetrosianI l|1976t) system (sec Blant on et all [2001 for a 
discussion of the advantages of Petrosian magnitudes), 
and are corrected for Galactic extinction usin g the dust 
maps of iSchlegel. Finkbeiner. fc DavisI l|1998t ). Further 
details about the photom etric data and the parani eters 
measured may be found in lLupton et all l)200lil2005|) and 
IStoue-hton et a,l. (2002V 

2.2. Cluster Finding Technique 

Clusters used in this study are detected by 
the maxBCG algorithm. This method relies on 
the observation that clusters host a population of 
early- type galaxies that have small dispersion in 
color o wer. Lucev. fc Elh sllT99l [Stanford et al.llT99l 
l^rail ct al.1 119981 IGladder s fc Yed 1200(1 . These clus- 
ter members populate the red sequence on a color- 
magnitude diagram. The brightest cluster galaxies 
(BCGs) have colors t hat are compatible with the 
red sequence gala xies fArago n-Salamanca et al.' '199^ 
INelson et al.ll200^ . and also have a very small disper- 
sion in lumin osity (Sandagc 1972, 1976; Oemlcr 1976; 
HS£22£^^^ ^80j ,yS chneider. Gunn fc Hocssel ISS^ 
Postman fc Laueilll995HCollins fc Ma!m^ il998'). There is 
evidence that the colors and magnitudes of BCGs remain 
predictabl e (and thus indicative of redshift ) until at least 
z = 0.6 ijAragon-Salamanca et alJ Il998t INelson et al.l 
12002(1 . shifting only due to passive evolution. 

The maxBCG algorithm takes advantage of these ob- 
served properties to find BCGs and the galaxies that 
are associated with them. For every galaxy, the algo- 
rithm calculates the highest likelihood, £max, that any 
given galaxy is a BCG (which means the galaxy has the 
properties of a BCG plus has a red sequence around it), 
and then identifies clusters by finding galaxies with high 
-Cmax values compared to the surroundings. >Cmax for 
each galaxy is determined by finding the maximum of 
C{z), the calculated likelihood as a function of redshift 
for the candidate. The likelihood function C{z) is the 
sum of two terms, and is calculated at every redshift 
from z — 0.0 to z = 0.6 in steps of 0.01. That is, for each 
galaxy we calculate 

-Cmax = max C{z); where £(z) = £bcg + logNgais. (1) 

The first term, £bcG7 is the BCG likelihood: the likeli- 
hood that the galaxy in question is consistent with, com- 
pared to the known population dispersions, the appar- 
ent magnitudes and colors of the mean BCG population 
as seen at that redshift. These colors and magnitudes 
were derived from the properties of Abell clusters ob- 
served by the SDSS. The second term, log Ngais, is the 
log of the count of the other galaxies in the vicinity that 
are also of the right color and magnitude to be cluster 
members. To count the number of galaxies, we examine 
only those projected within lh~^ Mpc (at the redshift 
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Fig. 1. — Color-magnitude diagram of observed g - r versus 
apparent j-band for galaxies near a rich cluster at 2 = 0.15. Ellipses 
represent 1-, 2-, and S-cr contours around the mean BCG color and 
magnitude at that redshift. The dotted line indicates the track of 
BCG color and magnitude as a function of redshift. The horizontal 
lines and vertical dashed line show the region of inclusion for Ng^ia 
determination. 



in question) that fall in the red sequence for each red- 
shift. A galaxy is within the red sequence if it is within 
2(7 (cr = 0.05 mags) of the mean BCG color at that red- 
shift, fainter than the BCG candidate, and brighter than 
Mi = -20.25 (approximately 0.5 L*). Figure ^ shows a 
color-magnitude diagram for a rich cluster. The region of 
inclusion for Ngais determination shown by the contours; 
the color and apparent luminosity of a passively-evolving 
BCG as a function of redshift is indicated by a dotted 
track. The interplay between finding a galaxy with the 
right color and luminosity to be a BCG and a galaxy 
with red sequence neighbors determines the redshift dis- 
tribution of C{z). We find the redshift that maximizes 
C{z) for the galaxy in question. That galaxy is assigned 
-Cmax and the corresponding redshift. At the end of this 
process each galaxy in the catalog has a single maximum 
likelihood and a single redshift. 

We then need to select the cluster centers out of the 
catalog of all galaxies. We find the peaks in the dis- 
tribution of £max over all galaxies; these peaks are the 
clusters. For each candidate BCG, we check whether 
its likelihood £max is the highest likelihood when com- 
pared with neighboring candidates within Az = 0.05 that 
are projected within lh~^ Mpc of the BCG candidate in 
question. If the candidate BCG's £max is the greatest 
of those neighbors, we list the candidate in the cluster 
catalog as a BCG. 

The cluster catalog produced contains information 
about each cluster identified, including the photometric 
properties of the BCG, the estimated redshift and the 
richness, Ngais- Each cluster center is taken to be at the 
location of its BCG. The richness measurement for the 
cluster, Ngais, is defined to be the number of galaxies in 
the red sequence at the derived redshift. The resulting 
catalog contains objects over a wide range of richnesses: 
from quite poor systems of only a few galaxies (Ngais < 
8; cr„ < 300 km s^^; 10,560 systems in 0.07 < z < 0.3) to 
very massive clusters of hundreds of galaxies (Ngais > 30; 
cr„ > 700 km s-^ 19 systems in 0.07 <z< 0.3). Figure 
[3 shows the distribution of identified objects as a func- 
tion of redshift and of richness. So as to avoid making an 
arbitrary distinction between a group and a cluster, we 
will generically refer to a system of galaxies as a cluster. 
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and specify the value of Ngais of that system. 

The maxBCG algorithm has been tested extensively 
for completeness and purity. All previously known Abell 
and NORAS X-ray clusters in the region surveyed are 
recovered. Simulations suggest that maxBCG recovers 
and correctly estimates the richness for greater than 
90% of clusters and groups present with Ngais > 15 
out to a redshift of z = 0.3. The completeness and 
selection function of the algorit hm will be further ex- 
plored in iWec hsler ct al.' f2005). The clusters identi- 
fied by this algorithm have been compared with the ob- 
jects found by different cluster- finding algorithms run 
on the same dataset. Discussion of the differences be- 
tween maxBCG and other algorithms can be found in 
iBahcall et aP f2003). 

One of the strengths of this algorithm is that it is 
a robust photometric redshift estimator for the clus- 
ters: for the 6708 clusters in the catalog with spec- 
troscopic redshifts available, the dispersion between the 
maxBCG estimated redshift and the spectroscopic red- 
shifts is Az = 0.018, as seen in Figure|31 and is smaller for 
the highest richness clusters. These 6708 systems span 
the range of Ngais in the catalog. As we do not have spec- 
troscopic redshifts for all clusters examined, but are con- 
fident in relying upon these estimates, we will henceforth 
use the term redshift to mean the estimated photometric 
redshift determined by maxBCG. 

3. BACKGROUND SUBTRACTION 

Given a set of cluster centers with well-defined three- 
dimensional positions, we need to find the galaxies asso- 
ciated with those clusters. This section explains how we 
apply background subtraction techniques to the SDSS 
data, and how we check our method by constructing and 
examining the radial density profile and the luminosity 
function. 

Generically, any properties of galaxies in clusters can 
be described by some population distribution function 
(PDF) in a multi-dimensional parameter space. Prop- 
erties of each galaxy such as luminosity, color, star for- 
mation rate, mass, and distance from the cluster center 
may be used as the parameters of a PDF. Examining 
how the galaxies of a particular cluster occupy the pa- 
rameter space is a way to sample the PDF; with a large 
enough set of cluster galaxies, we may statistically deter- 
mine the PDF quite well. The overall properties of the 
clusters (e.g. cluster mass or X-ray temperature) may be 
used to identify different sets of clusters, and the PDFs 
of galaxies in these different cluster samples compared. 
In this way, we can explore how the properties of galaxies 
in clusters are related to the characteristics of the host 
clusters. 

Without redshifts for all galaxies, we can only examine 
the PDF of galaxies associated with clusters by making 
an appropriate correction for the set of field galaxies pro- 
jected by chance along the line of sight to the clusters. 
We assume the presence of a cluster at some redshift 
does not affect field galaxies found along the same line 
of sight. That is, the PDF of galaxies projected around 
a cluster center has two independent components: the 
distribution of real cluster galaxies, and the distribution 
of random background and foreground galaxies. To de- 
termine the projected, azimuthally averaged PDF of just 
the galaxies associated with clusters (the PDFc), we ex- 



amine the PDF of all galaxies projected around cluster 
centers (the PDFcf) and the PDF of all galaxies pro- 
jected around a set of random (field) points on the sky 
(the PDFj-). The PDFc is determined by subtracting: 
PDFc = PDFci^ - PDFi?. Although we cannot iden- 
tify exactly which galaxies make up a particular cluster, 
we can very accurately describe the mean properties of 
galaxies associated with a set of clusters. 

There are a variety of ways in the literature for 
measuring the contribution of non-cluster members 
without having spectroscopic information. Historically, 
the population of field galaxies was esti mated from 
number-flux counts in sepa rate surveys fe.g. lAbeIilll958t 
iLuggedlToM lColles£llT98 9l). although this method has 
the disadvantage of not having the cluster and field sam- 
ples measured in t he same set of dat a. More recentl;^ 
some authors llValotto. Nicotra^ Muriel. & Lambaa 
19971 IPa ohUo e t all I2001| IGoto et al.. .2002 
Popesso. B oehring er fc VogesI 1200^1) have measured 
the background in an annulus centered on the cluster, 
in order to ensure that the background measurement 
is made using data of similar depth and seeing as the 
dat a in the region of the cluster. Other authors, such 
as lAndreon et al.l 1)2004(1 estimated the background 
from a nearby control field, or from the logN — logS 
relationship from the s ame dataset from whi ch the 
cluster sample is drawn llLin Mohr. fc StanfordI |2004) . 
IGarilli. Maccagni. &: AndreonI l)1999j) subtracted in- 
terlopers on the basis of color information, removing 
"galaxies with colors not matching the expected ones 
at the cluster redshift." As galaxy surveys increase in 
area, it becomes feasible to measure the background 
counts directly from the general field, as done by 
IGladders fc ^QQB)- The SDSS data offer large 
regions of sky measured to the same depth and with the 
same seeing, so we are able to determine the contribution 
of field galaxies in the same data as the clusters without 
artificially restricting the field measurement to the 
cluster neighborhood, or making assumptions about the 
color or luminosity distribution of cluster members. By 
using a set of random points as the locations around 
which the field galaxy population is determined, we 
measure the characteristics of all galaxies that are 
associated with clusters. 

3.1. Application to SDSS data 

The SDSS is an ideal dataset with which to examine 
the PDF of cluster galaxies because it provides sky cov- 
erage for a large number of clusters and ample blank 
sky for measuring the field distribution. The SDSS data 
offer a rich parameter space with which to define the 
PDF. Properties such as luminosity, color, star forma- 
tion rate, and morphology may all be explored. In this 
work, for g-, r-, i- and z-bands, we construct and exam- 
ine the PDFc (Ngais, M): the density function of clus- 
ter galaxies per surface area in a three-dimensional space 
of cluster richness Ngais, projected radius r, and absolute 
magnitude M. The PDF may then be projected onto the 
axis of absolute magnitude to show the luminosity func- 
tion of cluster members, or onto the axis of projected 
radius to show the radial density profile of the cluster. 
In this section we describe the samples of galaxies exam- 
ined, and as the PDFcf and PDF^ are constructed in 
the same manner, discuss the construction of a general 
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Fig. 2. — Left: the redshift distribution of maxBCG-identified objects. We use clusters in the shaded region; simulations show that the 
completeness rate begins to drop for z > 0.3. Right: the distribution of 0.07 <z< 0.3 maxBCG-identified objects as a function of cluster 
richness, Ng^is. There are 12,830 systems identified in this redshift range; 2270 of them have Ng^is > 8, and 19 clusters have Ng^is > 30. 




Spectroscopic redshift 



Fig. 3. — Photometric redshift estimation for the maxBCG clus- 
ter finding algorithm is tested here by comparison to spectroscopic 
redshift determination for a total of 6708 maxBCG clusters. The 
typical photoz error is cr^ = 0.02 for the full sample, falling to cr^ = 
0.014 for the redshift range 0.07 < z < 0.3. The increase in photoz 
errors around z = 0.37 occurs because this is where the 4000 A 
break, the most significant feature in a typical galaxy spectrum, 
passes from the g to r filters. The clusters used in the comparison 
span the full range of Ng^igOf the catalog. 



PDF(Ngais, r, M). 

3.1.1. Cluster and Field Samples 

To measure the PDFci?, we examine galaxies projected 
near the 12830 maxBCG objects found in the redshift 
range 0.07 < z < 0.3. These systems have richnesses in 
the range 2 < Ngais < 66. We use all galaxies projected 
within 2h~^ Mpc of the cluster centers and in the abso- 



lute magnitude range -24 < M < -16. We bin the data 
in 50 kpc radial x 1 Ngais richness x 0.2 absolute 
magnitude bins. The number of galaxies in each bin is 
then normalized by the physical area observed in each 
bin. Details of the calculations are discussed below. 

To determine the PDF/t-, we examine galaxies along 
lines of sight to randomly chosen field locations. For each 
cluster, we choose five positions on the same ~ 200 deg^ 
stripe of sky as the cluster, with random RA and DEC. 
These field positions are assigned the same redshift as 
the cluster, and labeled with the richness of that cluster. 
The resulting set of 64,150 field points are observed with 
the same seeing and to the same depth as the clusters, 
and for any set of clusters there is a set of field positions 
with the same redshift distribution. The same radial and 
magnitude ranges and bin widths used for the PDF cf 
are applied to determine the PDF^?. All excess galaxies 
seen around cluster locations as compared to the field are 
identified as cluster galaxies. 

3.1.2. Absolute Magnitudes 

To calculate absolute magnitudes M, apparent mag- 
nitudes m must be corrected for luminosity distance, 
Galactic dust extinction, and /^-corrections as 

M = m-5logio(^^^^^ -R-K{z), (2) 

where Dl{z) is the luminosity distance for our as- 
sumed cosmology ' R is the correction for reddeni n g, corn - 
puted following iSchleeel. Finkbeiner. fc Davis! l)1998j) : 
and K{z) is the appropriate if-correction. 

if-corrections are necessary because galaxy magni- 
tudes in the observed bandpasses correspond to different 
rest-frame magnitudes depending on the redshift of the 
galaxy. To compare magnitudes of galaxies at different 
redshifts, we apply a if-correction to convert all the mag- 
nitudes to a fixed set of bandpasses. To do so, we use 
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Fig. 4. — Left: Radial profiles of clusters for selected bins in richness. Shown are very poor (upper left), poor (upper right) medium 
(lower left), and high (lower right) richness clusters. The galaxies identified as BCGs have been removed. The lowest richness objects tend 
to be in under-dense regions. For clusters with Ng^ia > 8, a power law is a reasonable fit to the profiles. The best fitting power law is 
shown. Right: Best-fit power law parameters for the radial profiles of clusters of different richness. Shown are 1-, 2-, and 3-a contours 
of the model parameters for all richness bins with Ng^ig > 8. 



the meth od oflBlanto n et al.l l|2p03b: kcorrect v3_2f) . and 
following iBlanton^riiLrpOOSar . j^-correct all galaxies 
to z = 0.1. This redshift is chosen since it is close to the 
median redshift of the SDSS spectroscopic sample, and 
thus requires the smallest typical corrections. 

Although the true redshifts of the galaxies are un- 
known, we apply fC-corrections and calculate luminos- 
ity distances as though all galaxies projected around a 
given point (cluster center or field location) are at the 
same redshift as that point. For the galaxies distributed 
at different redshifts along the line of sight, the result- 
ing absolute magnitudes are not correct, but for galax- 
ies actually located at the redshift of the given position 
(i.e. those galaxies physically associated with the clus- 
ter) the X-correction and Dl{z) are appropriate. When 
the PDFiT' is subtracted from the PDFc_f, the contribu- 
tion from galaxies not at the redshift of the cluster is 
removed, leaving only the cluster galaxies for which Dl 
and K{z) are correctly determined. 

Our color-dependent 7^-corrections also affect the ab- 
solute magnitude to which the sample is complete and 
volume-limited at a given redshift. Galaxies of different 
colors (therefore with different amounts of if-correction 
applied) and different apparent magnitudes can have the 
same absolute magnitude. Thus if-correcting a range of 
uncorrected absolute magnitudes SMuncorr maps these 
magnitudes to the same Mkcorr- For example, at z = 0.3, 
6Muncorr ~ 0.4 mag; at that redshift a galaxy with ap- 
parent r-band magnitude of 21.0 (the survey limit) cor- 
responds, prior to iiT-correction, to an absolute magni- 
tude of -19.9. To be complete to z = 0.3, we can only 
use those Mkcorr for which the corresponding range of 
SMuncorr does uot cxtcud fainter than -19.9. The end 
result is that we adopt more conservative completeness 
limits to avoid color bias at faint luminosities. We do 
not use data in u-band because both i^T-corrections and 



star-galaxy separation are not as robust in this passband. 
The resulting absolute magnitude limits for a complete, 
volume-limited to z = 0.3 sample of galaxies in r, i, 
and z are thus -20.2, -19.6, -19.4, and -20.6 respectively. 

3.1.3. Effects of Geometry and Luminosity 

We correct for incompleteness both in geometry and in 
luminosity. Geometric incompleteness occurs when the 
search radius around clusters extends beyond the bound- 
aries of the survey. For example, at a redshift of z = 0.07 
(the lowest redshift cluster considered here), the 2h'~^ 
Mpc radius aperture is ~ 0.75°in diameter; some clusters 
he too close to the edge of our 2.5°-wide stripe of sky to 
have all galaxies within the desired aperture contained 
on an observed region of sky. For each radial bin of each 
cluster, we account for this geometrical incompleteness 
by calculating the area that lies on an observed region 
and weight the galaxy counts in each bin accordingly. 

Luminosity incompleteness arises because the apparent 
magnitude limit of the survey causes a varying range of 
absolute magnitudes to be accessible at varying redshifts. 
For example, to a redshift of 0.3, we can only see galax- 
ies with Mr < —19.6 but can examine all clusters in our 
catalog, while a,t z — 0.07, we can see to ~ —16.5, 
but are limited to only a few clusters. To avoid restrict- 
ing ourselves to studying only galaxies brighter than the 
completeness limit of the full set of clusters (to z = 0.3), 
we account for this luminosity incompleteness. For each 
magnitude bin, we determine the number of clusters at 
redshifts low enough to have galaxies observable to that 
limit. The galaxy counts in each bin are weighted accord- 
ingly. The result is that the bright end of the luminosity 
distribution is based on galaxies in all clusters in the 
catalog (at redshifts out to z = 0.3), but the faint end is 
determined from galaxies associated with lower-redshift 
clusters only. For our determination of i?200i 
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solute magnitude limits that ensure luminosity complete- 
ness for the full set of clusters, and push fainter only for 
examining the luminosity function of cluster galaxies. 

Having calculated the radial and absolute magnitude 
distributions of galaxies, we can determine the normal- 
ized PDF(Ngais,r,Af ) per surface area for galaxies around 
any set of positions and test our algorithm. 

3.2. Consistency Checks 

In this section we check the background subtraction 
technique by examining the radial distribution and the 
luminosity distribution of galaxies around field loca- 
tions, and compare these distributions to those measured 
around cluster centers. 

Note that only a background-subtracted PDF (e.g. the 
PDFc) contains physically meaningful information; the 
PDF around any set of points before subtraction is dom- 
inated by galaxies projected by chance along the same 
line of sight, which have not been properly /C-corrected. 
We therefore reserve the name 'luminosity function' for 
the projection of a background-subtracted PDF onto the 
axis of absolute magnitude; such a projection for a non- 
subtracted PDF we refer to as a 'luminosity distribu- 
tion.' Likewise, we will reserve the term 'radial pro- 
file' for background-subtracted PDFs only, and will refer 
to the 'radial distribution' of galaxies when discussing a 
non-subtracted PDF. 

A check on the errors recovered by our background 
subtraction technique is done by comparing the PDF^? 
with the population density of galaxies around a set of 
random points, which are different random points than 
those used in constructing the PDF^?. The PDFs mea- 
sured around two different sets of random points should 
be statistically identical. 

3.2.1. Radial Number Density Profile 

A radial number density distribution is constructed by 
projecting a PDF(Ngais, r, M) onto the axis of projected 
radius. Examining the distribution of galaxies around 
field locations enables us to check our correction for in- 
completeness due to geometry. 

When we construct the radial profile of clusters, we 
first examine the profile using only galaxies brighter than 
the completeness limit for the full set of clusters, and 
then using all galaxies while accounting for the varying 
completeness limit as described above. The results are 
statistically the same, but allow for inclusion of fainter 
galaxies in the latter case. 

The radial distribution of galaxies around random 
points is flat, demonstrating that we are properly cal- 
culating the area observed when correcting for geomet- 
rical incompleteness. Comparing the radial distribution 
around field locations with that measured around a dif- 
ferent set of random points, we see that the distributions 
are statistically identical at all radii. Any differences are 
within the error bars, which reflect the Poisson fluctua- 
tion of our sample. 

Around cluster centers, however, we find a significant 
excess of galaxies compared to the field. This excess 
varies as a function of radius and as a function of richness. 
The left-hand panel of Figure 0| shows the radial profile 
for selected richness bins corresponding to very low, low, 
medium, and high richness clusters, with the galaxies 
identified as BCGs removed. For clusters with Ngais > 5, 




Fig. 5. — Luminosity function per unit surface area in r-band for 
galaxies in rich clusters (30 < Ngais < 66), using galaxies within 
a fixed 1.5/i~^ Mpc aperture. Top (before subtraction): lumi- 
nosity distribution of galaxies projected around cluster centers and 
field points (dotted line). Because we /f-correct all galaxies to the 
cluster redshift regardless of the true redshift of each galaxy, the lu- 
minosity distribution around random points does not look l ike the 
luminosity function of galaxies as measured by e.g. Bl anton et alj 
1 2003a). Bottom (after subtraction): luminosity function of 
galaxies associated with these rich clusters. The solid line is the 
best-fitting Schechter function. 



a power law is an acceptable fit to the profiles. For sets 
of clusters of different richness and Ngais > 8, the radial 
profile is roughly consistent with ^ 7-^11 surface density 
profile, and thus a r^^'^ volume density profile, but with 
increasing normalizations for richer clusters, refiecting 
the larger size of more massive clusters. For all richness 
bins above Ngais > 8, the contours of the power law 
parameters are plotted as a function of cluster richness 
in the right-hand panel of Figure 0] In [jS] we further 
examine the radial profiles, using our measurement of 
i?^o to fit aiNavarro, Frcnk. & White (1997) profile. 

The lowest richness objects (2 < Ngais < 4) tend to be 
in under-dense regions. That is, compared to the distri- 
bution of galaxies around a random point in the universe, 
these systems of just a few red galaxies tend to have few 
nearby neighbors. We do not expect these very low Ngais 
objects to be representative of the full population of low 
mass halos, as selection effects of the niaxBCG algorithm 
are significant for these systems. Nonetheless, such a 



8 



Hansen et al. 




-24 -23 -22 -21 -20 -19 -18 -21.2 -20.8 -20.4 -20.0 

M,. M. 

Fig. 6. — Left: Luminosity function per unit surface area in r-band for clusters with 18 < Ng^ig < 66 in three redshift slices. Cluster 
members used are projected within l.5h~^ Mpc of the cluster center. The BCGs have been removed. Squares are for clusters in 0.07 < z 
< 0.15; diamonds for 0.15 < z < 0.20; stars for 0.20 < z < 0.30. The best-fitting Schechter functions are overlaid. Right: The 1-, 2-, and 
3-0" contours of the Schechter function parameters. The shift in the characteristic luminosity, M, , reflects the passive evolution of the 
early-type galaxies. 



sample is interesting; we discuss very low Ngais systems 
in more detail in S0 

3.2.2. Luminosity Function 

The luminosity function of galaxies in clusters is the 
projection of the PDFc(Ngais, r, M) onto the absolute 
magnitude axis. The result is the mean number of galax- 
ies per cluster per unit surface area as a function of lu- 
minosity and of richness. 

We check our background subtraction and incomplete- 
ness correction by comparing the luminosity distribution 
around two different sets of random points, expecting 
no significant difference. We begin by selecting only 
locations with low redshift, so that the luminosity dis- 
tribution may be examined to faint magnitudes with- 
out weighting for the redshift distribution, and then also 
check the distribution with the full sample and appro- 
priate weighting as discussed in The two samples 
have identical distributions. 

The luminosity distribution contains a statistically sig- 
nificant excess of galaxies around cluster centers com- 
pared to galaxies around field locations. Figure |5l shows 
our determination of the luminosity function of galax- 
ies that are in rich clusters (30 < Ngais < 66) and are 
within 1.5h~^ Mpc of the cluster center, both before and 
after subtraction of the field (top and bottom figures re- 
spectively). In the bottom figur e the solid line p lots the 
best-fitting Schechter function (!Schechte3ll97fiD . of the 
form 

(f>{M)dM = 

0.41n(10)(/.aO-0-4(*'^-*^*)("+i)e-i"""*"""*'dM (3) 

where a is the faint-end slope and A/* is the turnover 
magnitude; we fit the data using the Levenberg- 
Marquardt minimization procedure. We recover a 
luminosity functi on that is compar able to that of rich 
clusters found bv lGoto et al.l l)2002|) . who used different 



cluster-finding and background subtraction algorithms 
with the same sample of SDSS data used in this paper. 

We also examine the LF of cluster galaxies in three 
redshift slices to test whether our weighting scheme for 
the faint end is correct. The left panel of Figure El shows 
the LF for clusters with 18 < Ngais < 66 for 0.07 < z < 
0.15, 0.15 < z < 0.20, and 0.20 < z < 0.30. Galax- 
ies projected within 1.5h~^ Mpc are used; the BCGs are 
not included. We fit each distribution with a Schechter 
function (solid lines). The confidence ellipses for the fit 
parameters a and are plotted in the right panel of 
the figure. The only difference we detect between the 
LFs in different z slices is at the bright end, as reflected 
by the shift of toward fainter magnitudes at lower 
redshifts. This shift is comparable to what is expected 
due to passive evolution of the red sequence galaxies in 
the clusters. As our if-corrections are non-evolving, we 
expect A'U to be 0.25mag brighter for the highest red- 
shift bin than for the low z bin due to passive evolution. 
We note that these LFs are measured within a fixed phys- 
ical aperture for systems of a wide range of richnesses. 
We present this comparison of luminosity functions for 
redshift slices only as a check that we are recovering sen- 
sible LFs. A more detailed investigation of the evolution 
of the LF of cluster galaxies will be done in later work. 

Since more massive clusters are larger, measuring the 
LF (or any other projection of the PDF that varies ra- 
dially) within a fixed physical aperture samples different 
parts of clusters of different richnesses. For example, as 
can be seen from the radial profiles in Figure^ a lh~^ 
Mpc radius around a poor group encompasses the en- 
tirety of the group, but only samples the inner region of 
a rich cluster. Thus, in order to compare the LF of clus- 
ter galaxies in clusters of different richness, we should 
examine the LF of only those galaxies within some aper- 
ture that scales with richness. In addition, the luminos- 
ity function of cluster galaxies may vary with radius. To 
compare radial trends in clusters of different richness, we 



R200, Radial Profiles and Luminosity Functions for SDSS Galaxy Clusters 



9 




TABLE 1 

Power Law Fits for -R2oo{^gais) 



80 100 



Fig. 7. — Relationship of the cluster characteristic radius, iJ2oo' 
here measured in r-band, to richness, Ng^is- Diamonds mark the 
data using the over-density threshhold AMcrit = 200 times the 
critical density; the other data points use the over-density thresh- 
hold AA/" mean = 200 times the mean density. Using either over- 
density threshhold results in a trend well fit by a power-law with 
the same scaling. TablelTllists the best-fit values of the parameters. 



should also use an appropriately scaled aperture. In 21 
we present our calculation of a characteristic radius of 
clusters as a function of richness, and will return to a 
discussion of the dependence of the luminosity function 
on richness and radius in SJS] 

4. iJ^o DETERMINATION 

To appropriately compare properties of low- and high- 
mass objects we need to understand how the characteris- 
tic size of clusters varies with richness. Motivated by the 
way i?A is defined in N-body simulations, as the thresh- 
old radius interior to which the mean mass density of a 
cluster is A times the average mass density, we define 
an analogous R%q using the space density J\f of galax- 
ies in clusters compared to the average space density of 
galaxies. If galaxies are unbiased with respect to dark 
matter on all scales, we would have = Ra- Since 
the bias is close to unity, we accept i?^ as a reasonable 
approximation. Following simulations, we use A = 200 
as our threshold mean over-density of cluster galaxies, 
which occurs at the radius i?2do- 

Some authors take the average mass density to be the 
critical density, while others use the actual mean back- 
ground density. For the main result of this paper, the 
scaling of i?2oo with cluster richness, we present results 
both using an over-density threshold of AAfcrit = 200 
times the critical density and using AAfmean = 200 times 
the mean background density. We find the same scaling 
using cither threshold. Wc intend to compare this work 
with the resu l ts usin g the Hubble Volume simulations of 
lEvrard eTal] pOO^ . who use an over-density threshold 
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measured with respect to the critical density. Therefore, 
for investigations regarding the cluster galaxy population 
within R^Q, we present results using AMcrit- Through- 
out this work, we use the term i?2oo to mean the radius 
interior to which the mean number density of galaxies 
is 200il^ times the mean space density of galaxies, or 
equivalently, 200 times the critical density. 

To determine the mean space density of field galax- 
ies, we use the r, i, and z luminosity functions of 
iBlanton et al. (2003^, which are properly normalized to 
a volume density and are determined with an SDSS spec- 
troscopic sample of galaxies from the same region of sky. 
We integrate these field LFs down to the absolute mag- 
nitude hmits apphed to our cluster sample (-20.2, -19.6, 
-19.4, and -20.6 in g, r, i, and z respectively), and take 
the resulting value to be the average space density in 
that passband. We use only these four bands since the 
u-band if-corrections and star-galaxy separation are not 
as robust as in these bands. 

To measure the mean space density of cluster mem- 
bers, we use the PDF(7(Ngais, r, M) determined above 
for each richness. For clusters of a given richness and in a 
given bandpass, in each radial bin r we sum over all bins 
with radius < r and with absolute magnitude brighter 
than the completeness limit. We assume the galaxies are 
contained in a sphere of radius r to calculate the volume 
density, then divide by the mean space density of field 
galaxies to get the fractional excess. In actuality, the 
galaxies are contained in a cylinder of diameter 2r. Thus, 
although the uncorrelated galaxies are removed from the 
measurement, there is an excess of galaxies at radius r 
due to the projection. This excess depends on the shape 
of the radial density profile. For samples that have the 
same radial profile, this excess is simply the same multi- 
plicative factor for all. For example, if the profile is l/r^, 
the factor is tt/2. 

We then use the binned mean over-density vs radius 
information to find the radius interior to which the num- 
ber density of cluster galaxies is 20057"^ times greater 
than the field density. Since the density has roughly a 
power law radial dependence, to determine exactly where 
ATV = 200ri~^, we fit a fine in log- log space in the region 
in which the over-density passes through 20017"^. The 
errors on this fit include the uncertainty on the density 
value and the width of the radial bin, and determine the 
uncertainty in the i?200 value. In this manner we calcu- 
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late R2Q0 for each Ngais bin for which there is at least one 



The relationship between characteristic radius i?200 
and cluster richness Ncrais is well fit by a power law with 



10 



Hansen et al. 




Fig. 8. — Number density profiles for satellite galaxies in clusters of different richness. The best-fit projected NFW profile is shown, 
obtained by assuming the previously measured relation between i?2oo ^'^'^ Ngais, and fitting the data within i?200 °'^'y- T^^^ ^200 
responding to the log mean Ng^ia value of the richness bin is listed. The richness bins are the same used in examing the radial profiles 
previously {Figure|4j. The excess beyond i?2oo t° expected as the two-halo term becomes important. 



index ^0.6. The determination of i?200 ^^r clusters as 
a function of richness is the principle result of this work. 
Figure [3 shows R^q measured with r-band dat a as a 
function of Ngais (diamonds), with the best fit power law 
plotted. We also plot the relationship between radius 
and richness measured using AAfmean- The scaling is 
the same, with different normalization. Similar results 
are obtained for g, i, and z; the best fit power law pa- 
rameters with 1-(T uncertainties are listed in Tabled for 
all four passbands and both over-density thresholds. 

Under the assumption that the i?^g we have measured 
here for clusters is a good proxy for a mass density based 
R200 for dark matter halos, what sort of scaling relation 
would we expect with Ngais? A detailed answer to this 
question requires understanding the mass-to-light ratio 
as a function of both cluster mass and cluster radius, 
which is beyond the scope of this paper. In the simple 
scaling arguments below, we assume that these are both 

constant. The radius and mass of a cluster scale as i?2oo 
1/3 

~ M20Q , and the number of galaxies within i?2oo is likely 
to scale as a power law with cluster mass as NgaZ.flaoo ~ 
M"nn - This power h as been found to be close to unity 
(e.g., Kravtsovet al. 2004: Lin. Mohr. fc Stanfordll2004l: 
IZehavi ct al. 2004; Wcchsler ct al. 2005). The Ngais we 
use here, however, is measured within a fixed lh~^ Mpc 
radius aperture; it will typically be smaller than Ngai_R2oo 
for the most massive clusters and larger than Ngai_R2oo 
for smaller groups and clusters, and roughly consistent 

with Ngais ^gai.fiaoo' ^^^^ (3 < 1. If the galaxies fol- 



low an NFW profile out to max(l Mpc, i?20o), with a 
concentration around 5, one would expect (3 to be in the 
range ^ 0.50 — 0.65. For maxBCG clusters that have 
been found in the simulations of Wcchsler ct al. (2005), 
we find something similar: Ngais ^ ^ga/^flaoo' P^^ffii^g 
this all together, we have 

^200 - M200 - ^gal_R2oo ^gals ^ ^ gals, W 

which is in excellent agreement with the scaling relation- 
ship that we find for i?200 N^gais- This comparison 
suggests that our observationally determined i?200 a 
reasonably good proxy for i?200- Note that in detail the 
relation between Ngais measured at a fixed radius and 
Nga/_R2oo is not expected to be a power law over all halo 
masses, which implies that the power law relation found 
between i?200 and Ngais may break down when measured 
over a wide range of halo mass. This relationship, as 
well as the detailed relationship between i?2oo and i?2oo 
as traced by dark matter, will be explored further in 
future work, using a larger sample and comparison to 
simulations. 

5. GALAXY DENSITY PROFILES WITHIN 

Using our empirically measured i?2ooi we now exam- 
ine the radial density profiles of galaxies in these clusters 
in greater detail. Simulations suggest that dark matter 
halos have mass profiles characterized by a scale radius 
Vs = i?20o/cDj\/, where cdm is the concentration parame- 
ter for the dark matter (NFW; iNavarro. Frenk. fc White! 
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11997(1 . Galaxies do not necessarily trace the same de- 
tailed distribution as the dark matter. In particular, 
for any given sample of galaxies chosen with some selec- 
tion criteria, a range of processes (e.g., dynamical fric- 
tion, tidal stripping, enhanced or suppressed star for- 
mation) may affect the distribution of those galaxies 
within their host dark matter halos. Still, several re- 
cent studies have suggested that the number density 
profile o f galaxies is well described by the NFW func- 
tion (e.g. Carlberg et al.'1997'; van der Marel ct al. 20001 
[Mahdavi fc GcUer 2004; Katgert, Biviano, & Mazurl 
l2004l:"lLin. Mohr. fc Stanford 2004); here we fit our radial 
profiles with the projected NFW profile and examine its 
dependence on cluster richness. 

We express the number density profile in three dimen- 
sions as n(x) = nQX~^{l + x)~'^, with normalization rig 
and X =Cgai t/R^^^q, where Cgai is t he concentration pa - 
rameter of the galaxies. Following i BartelmannI l)1996fl . 
we write the projected surface density NFW as the inte- 
gral 

S(x) = cos0(cos0 -I- '^r^de. (5) 

Cgal Jo ^200 

For several bins in richness, we express the radial pro- 
file in units oir/R^Q, where each cluster has been scaled 
by the i?2oo appropriate for its Ngais value, as measured 
in the previous section and specified in Tabled We fit 
the resulting number density profile within i?2oo with the 
profile specified in Eq. the results are shown in Fig- 
ure IHl The fit does well for most of the richness bins. 
There is some excess outside the virial radii, which is to 
be expected as the two-halo term begins to contribute to 
the distribution. However, for the low-richness systems, 
some of this excess may be due to a misidentification of 
the cluster center. We discuss this issue further in 

The variation of the concentration parameter with 
richness is shown in Figure El Note that we have de- 
fined the concentration parameter for galaxies with re- 
spect to the measured -R200 of the galaxy profile; if this 
is not equivalent to the dark matter -R200 then Cgai will 
change accordingly. In the top of Figure [HI we show the 
measured Cgai divided by the expec ted dark matt e r con- 
centration Cdm, using the model of lBullock et alJ 1(200 If) 
assuming a cosmology with Vim =0.3. and cts = 0.9. In 
order to make the comparison, we assume that R^q = 
Riooy which may not be the case, especially for distribu- 
tions with such different concentrations. In particular, 
preliminary indications from both simulations and weak 
lensing measurements indicate that R'2qo< ^200- Because 
the theoretical prediction for cdm{M) gets steeper at 
higher masses, in this case Cgai{M) /cdm{M) would be 
closer to a constant with mass M. 

It is clear that the profiles of galaxies in maxBCG 
clusters have significantly lower concentrations than the 
dark matter profiles measured in CDM N-body sim- 
ulations. This finding is in agreement with previous 
work Ce.g. iLin. Mohr. fc Stanford 207)1: iCarlberg et all 
119971 Ivan der Ma,rel et a1 .1120001) . which has found of Cgai 
~ 2 — 4 for cluster galaxies. These low values, how- 
ever, should be interpreted as indications of how galax- 
ies are distributed within dark matter halos, and not 
of the concentration of the dark matter of these ha- 
los. This point was emphasized by iNagai fc Kravtsovl 
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Fig. 9. — Best-fit concentration parameter of the galaxy profile, 
assuming that -R200 given by the previously determined relation 
between -R200 ^^'^ Ng^is. The error bars are determined from the 
1-(T region of the surface. The top axis shows AI200, assuming 
the same -R200 to Ng^is scaling, and under the assumption that 
-R200 ~ ^200 and rirn = 0.3. 



12.0041) in their investigation of a set of hydrodynamic 
simulations of clusters, in which they found values of 
Cgal ~ 2 — 7 for eight clusters where the dark mat- 
ter concentrations were 6 — 16. In general, for any 
population of galaxies in a host halo, the radial distri- 
bution is dependent on the dynamical and star forma- 
tion histories of the galaxies once they enter the host 
halo, and may depend sensitively on ho w the popula- 
tion is selected JPiemand et al. 2004; Gao et alJ 120041: 
INagai fc Kravtsovl l2004|) . iMandelbaum et a'Lri|2004D 
have found that low values of Cgai quite similar to what 
we measure here are required to match the galaxy-galaxy 
lensing observations in SDSS. 

6. CLUSTER LUMINOSITY FUNCTIONS FOR Mr < -18 
WITHIN 

We can use i?^Q to compare commensurate regions 
within clusters of different richness. We first measure 
the luminosity function per unit surface area of cluster 
members within the appropriate i?200 for each richness, 
then rescale by the area contained within i?200 to deter- 
mine the actual number of galaxies of each brightness 
within i?^o for each richness. We combine the LFs into 
six bins of Ngais. The top panel of Figure [TUl shows the 
LFs of all galaxies within i?2oo ! the bottom panel those 
galaxies identified as BCGs have been removed. 

With the BCGs removed, a Schechter function pro- 
vides a reason able fit to the data in all Ngais ranges. 
iBlanton et all lj2004 'l have examined the luminosity func- 
tion of low luminosity galaxies, finding that surface 
brightness selection effects bias the LF to lower values 
fainter than Mr ~ -18. That is, low luminosity, low sur- 
face brightnes s galaxies tend to be missed. In addition, 
IBlanton et alJ 12004 ) find that the field galaxy LF turns 
up fainter than Mr ^ -18, and that the shape of the LF 
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26 < N.^, <66 
18 < N°^, <25 
14 < N,., <17 
11 < N^^, <13 
8 < N„,. < 10 




TABLE 2 

SCHECHTER FiTS FOR r-BAND LUMINOSITY FUNCTIONS WiTHIN -R2OO' 

BCGs Excluded 



-24 -22 -20 -18 

Fig. 10. — Luminosity function within iJ2ooi expressed as the 
number of galaxies per magnitude, for galaxies in clusters of differ- 
ent richnesses. Top: LFs of all galaxies within R^q- A Schechter 
function is not a good fit to many of the richness bins, as the effect 
of the BCG is significant. Bottom: LFs for the same sets of clus- 
ters as above, but with those galaxies identified as BCGs removed. 
Schechter functions fit for Air < — 18 are shown; the faint end slope 
was held fixed for bins with Ng^is > 8. The parameters and oi 
the fits are listed in Table |21 



over a wide range of magnitudes is best fit by a dou- 
ble Schechter function. In hght of the concerns about 
missing low surface brightness galaxies at faint magni- 
tudes, we restrict our fit to Mr < -18. For Ngais> 8, 
a = —1 provides a reasonable fit to the data in this mag- 
nitude range. The parameters of the Schechter function 
fits, with a fixed, are listed in Table El The two lowest- 
richness bins are not well fit with a faint-end slope of -1; 
in the Table we list the best-fitting Schechter function 
parameters for these two bins. Detailed investigation of 
the differences in LF between different richness samples 
will be done in future work with a more extensive data 
set. However, we see that the primary change is the in- 
crease in normalization of the LF with increasing cluster 
richness, that the characteristic luminosity, Af,, bright- 
ens moderately toward richer clusters, and that the Ngals 
< 8 systems are different from the richer objects. The 



Richness 



0. 



xVd.o.f. 



2 < Ng,i,< 4 
5 < Ng,i,< 7 
8 < Ng3,i,< 10 
11 < Ng,i,< 13 
14 < Ng,i,< 17 
18 < Ng,i,< 25 
26 < Ngai,< 66 



-19.95±0.12 
-20.06±0.08 
-20.65±0.04 
-20.68±0.04 
-20.58±0.05 
-20.70±0.04 
-20.86±0.05 



-0.38±0.12 
-0.55±0.07 
-1.00 (fixed) 
-1.00 (fixed) 
-1.00 (fixed) 
-1.00 (fixed) 
-1.00 (fixed) 



0.94±0.07 
1.89±0.12 
2.93±0.09 
5.29±0.17 
9.53±0.33 
14.1±0.47 
29.9±1.07 



1.2 
1.2 
1.8 
2.3 
2.0 
1.1 
1.9 



mean redshift of the different richness cluster samples 
changes by A z ~ 0.03, so we do not expect significant 
luminosity evolution (AM, < O.lmag). We note that 
since these LFs are determined using cluster members 
that reside within R^q, the local density is the same on 
average for these galaxies, so the richness dependence of 
the LFs may not be attributed solely to variations in the 
local environment. 

The sensitivity of the shape of the LF to the cluster 
richness makes comparison between different catalogs of 
clusters difficult. Different definitions of richness and/or 
different bins of richness will result in different measured 
LFs. In addition, other catalogs typically present results 
for the luminosity function using galaxies within a fixed 
physical aperture, rather than within an aperture that 
scales with mass. Some authors rescale individual clus- 
ter LFs by cluster richness before creating a composite 
LF, but still examine the LF within a fixed physical aper- 
ture. Nonetheless, we do find qualitatively similar results 
to other authors. The luminosity function we find for rich 
clusters is similar to the LF of rich clusters presented in 
iCxoto et alJ l|2002D . who used the same SDSS data, but a 
different cluster-finding algorithm and different method 
of background subt raction. Our results are also in agree- 
ment with those of lPopesso. Boehringer fc VogesI l|2005f) 
(PVB). Like PBV, we see that the faint end of the LF 
picks up below Mr ~ -18, even though we are likely miss- 
ing some of these faint galaxies due to surface brightness 
selection effects. Our LFs in other bands are also com- 
parable to the measurements of other authors, who typi- 
cally find a steeper faint end in bluer bandpasses. Table 
1 of PBV lists the Schechter parameters for composite 
cluster LFs retrieved from the literature for a variety 
of bandpasses. The LFs in lower richness bins are also 
comparable with what has been found by other authors. 
The 8 < Ngals < 10 groups are comparable to those of 
iMartfnez et al.l H200 2) . with velocity dispersions ~ 300 
km s""'^; we find a similar result to theirs for the lumi- 
nosity function of these groups when we, like they, in- 
clude the BCG. The LFs of very low richness systems 
(Ngals < 8) have a falling faint end slope, and a bright 
end that is dominated by the BCGs of these objects. 
These are the same systems that are preferentially found 
in under-dense regions, as shown in ^3.2.11 We discuss 
these systems further in ^ 

With the BCGs included, a Schechter function is not 
a good description of the data, except for in the very 
richest clusters where they are only a small contribution: 
the BCG population adds a bright-end "bump" to the 
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LF of the other cluster members. The BCGs tend to 
be increasingly bright in higher-mass clusters, but also 
become less important to the total light with increas- 
ing cluster mass. We fit the LFs within iJ^g, BCGs 
included, with a Gaussian for the BCGs plus a Schechter 
function for the non-BCG galaxies. This model is a good 
fit to the data for all Ngais ranges. The top panel of Fig- 
ure ^2 shows the relative amplitudes of the Schechter 
function and the Gaussian function, fsch,tn evaluated at 
the mean of the Gaussian component as a function of 
mean log Ngais. Over this range the trend is linear, with 
fsch,t2 ~ Ngais/38.3. This scaling is the observed analog 
for the SDSS t o the conditional ba ryonic mass function 
investigated bv lZheng et al.l l)2004|) (see their Figure 9). 

We plot the mean luminosity of the Gaussian as a func- 
tion of mean log Ngais in the bottom panel of the Figure. 
The trend of brighter BCGs in richer clusters is evident. 
Over the range we probe, our data are consistent with 
a single power law, Lc ^ Ngais°'^. This result is similar 
to, but slightly steeper than, the scaling with mass found 
by other authors LLin fc M ohr 2004; Zheng et al. 2004; 
Ivan den Bosch et alJl2004j) . but note that we are plotting 
the scaling as a function of Ngais and not mass. Our data 
are all for multiple-galaxy systems, and as such cannot 
constrain the very low mass end of the distribution stud- 
ied by others. However, in order to facilitate compari- 
son with previ ous work l|Va le & Ostrikcr 2004, and the 
compilation of iCoorav fc Milosavlievic. 1,200^) ). we fit a 
double power law function with several parameters fixed 
to agr ee with the low-mass behavior of I Vale & Ostrikeil 
(^003). We find that the relationship between mean BCG 
luminosity and Ngais is consistent with a model of the 
form 



{NgaljNcf 
[l+{Ngals/Nc)V' 



(6) 



where the mean BCG luminosity, L^, is scal ed to L^, of 
the lum inosity function for SDSS galaxies (B lanton et al.l 
I2003a|) . The best fit parameters are Nc = 0.98 ±0.04 and 
7 = 3.52 ± 0.02, but we stress that these are degener- 
ate with the parameters we have held fixed; our current 
data cannot put a strong constraint on the exponent in 
the numerator (specifying the scaling for low masses). 
We note that the observation that BCGs are drawn 
from a different luminosity function from other cluster 
members is well known (see, for example, the review by 
iCollins et~al.<.20 03 and references therein). The trend of 
increasing central galaxy luminosity with cluster richness 
has also been noted in many other observational stud- 
ies, such as Sandage & Hardv (1973); Sandagc (1976); 
iHocssel ct al. (1980); S chneider. Gunn & Hoessell (|1983f) 
and Lin fc Mohr (20041) . and is consist ent with pre- 
diction s fr om the semi-analyt ic models of iBenson et al.l 
ll2rif)l and lZheng etlTI ll200l (Z04). 

Local density is known to correl ate with several galax y 
properties, including luminosity ijBlanton et al.ll2003aj) . 
and so we also expect to see differences in the LF as a 
function of radius. Figure El shows the LF of galaxies 
in three radial bins: 0.0 < r/R^Q < 0.25 (light lines), 
0.25 < r/i?^o < 0.75 (medium lines), and 0.75 < r/i?^^ 
< 1.75 (heavy lines) for galaxies in clusters in four bins 
of richness. For the innermost radial bin, we plot the 
LF both with and without (dotted lines) the BCGs. For 
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Fig. 11. — Parameters related to fitting the luminosity functions 
within i?200' BCGs included, with a model comprised of a Gaus- 
sian (for the BCGs) plus a Schechter function (for the non-BCGs), 
plotted as a function of the mean log Ngais for the ranges used. 
The LFs used are the same as in Figure ITOl top panel. Top: the 
ratio of the amplitudes of the Schechter and Gaussian functions 
evaluted at the mean fi of the Gaussian component. The best-fit 



line is plotted, and has fg 



ch,fi ' 



N 



gals 



/38.3. Bottom: The mean 



luminosity of the Gaussian fit to the LF of the BCGs, in terms of 
L* for SDSS galaxies, as a function of cluster richness. The solid 
line is L = Ngals^'''^; the dashed line is the double power law model 
of Equation IdI 



these LFs, we plot the luminosity function per unit sur- 
face area to explicitly show the change in local density. 
For all richnesses, the overall normalization of the LF 
decreases toward larger radii, as the density of galax- 
ies drops. Very bright galaxies are found predominantly 
only in the centers of both poor and rich objects. Ex- 
cept for the lowest richness systems, the faint end slope 
is roughly similar for all radial bins, although shifts 
somewhat fainter towards the outskirts of the clusters. 
The data are not statistically powerful enough to put 
strong contraints on these variations; future work with 
a more extensive data set will allow more detailed in- 
vestigation. We have not plotted the data for Ngais < 5 
objects, as just beyond i?2oo of these objects the radial 
profile becomes negative, as these objects live in under- 
dense regions. However, the LFs of these systems, both 
with and without BCGs, within 0.25r/i?2Qo is similar to 
the LF found for 5 < Ngais < 7 groups, although with a 
slightly fainter centroid of the BCG population. 

Our cluster-finding algorithm identifies the BCG as 
the galaxy that maximizes C{z)^ a statistic which in- 
corporates luminosity, color, and the number of fainter 
neighbors of similar color. There is an interplay between 
maximum local density and finding a galaxy with the 
right luminosity and color to be a BCG. Sometimes the 
maximum total likelihood selects a galaxy as BCG that 
the eye would not. The cluster center is defined to be at 
the location of the algorithm-defined BCG, and thus may 
not be at the true center of the cluster potential if either 
the most massive cluster member is not the one with the 
highest BCG likelihood, and/or is not at the center of 
the potential. Tests with simulations suggest that the 
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Fig. 12. — Luminosity function per unit surface area of galaxies in clusters of different richness as a function of distance from the cluster 
center in terms of i?200' Light lines are galaxies within i'/R'2qq < 0.25 (the dotted line indicates the LF with the BCGs removed); medium- 
weight lines are galaxies in the range 0.25 < r/R^Q < 0.75; the heavy lines are galaxies in 0.75 < i'/R'2oo < 1.75. For all richnesses, the 
overall normalization of the LF decreases toward larger radii, as the density of galaxies drops. Very bright galaxies are found predominantly 
only in the centers of both poor and rich objects. Note that the BCG population tends to be brighter and less dominant as cluster richness 
increases. 



centering of the maxBCG algorithm is good within ^ 
80ft.~^ kpc. For a rich cluster, this amount is a small 
fraction of the virial radius, but for a poor group, the al- 
gorithm may misplace the cluster center by a third of the 
virial radius. Even if the galaxy identified as the BCG is 
both the brightest and the closest to the cluster center. 
X-ray observations suggest that a BCG can reside as far 
as ~ 70/t~^ kpc from the center of the X-ray emission 
ijLazzati fc Chincarinil [T998: Lin fc Mohill2004|) . 

The inner region of the lowest richness objects may 
be biased due to failing to correctly identify the galaxy 
most closely located at the center of the cluster, and 
therefore incorrectly positioning the center of the clus- 
ter. For the radial profiles, a higher inner radial bin due 
to incorrect removal of the BCG of the lower-richness 
objects will skew the measured concentration parameter 
to higher values for those objects. We also see the pres- 
ence of BCG-like galaxies in the very inner regions of the 
lowest-richness clusters in the luminosity function: ex- 
amining the region within 50/i~^ kpc of the centers of 



poor clusters shows a luminosity function with a small 
"BCG bump" at the bright end, even when the maxBCG- 
identified BCGs are removed. Further investigation is 
needed to show whether these galaxies are the true BCGs 
of these groups, or whether the groups host a population 
of BCG-like galaxies in addition to the BCG. 

7. DISCUSSION AND CONCLUSIONS 

By using statistical background subtraction, we have 
observationally determined the changes in the radial dis- 
tribution of cluster galaxies as a function of cluster rich- 
ness, and used that information to calculate a charac- 
teristic radius for clusters of each richness. This model- 
independent radius, i?2ooi is based on the number density 
of galaxies analogously to the way i?200 in use in simu- 
lations is based on mass density. We find that i?^o 
hibits a power law scaling relationship with cluster rich- 
ness. This result is in good agreement with the expected 
scaling from simulations. 

Our result is also in general agreement with results ob- 
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tained by other groups usin g different methods fo r mea- 
suring richness and i?200- iCirardi et alJ l)199fiD exam- 
ined spectroscopically confirmed members of 90 rich clus- 
ters, and determined virial and core radii by fitting the 
observed galaxy distribution with an isothermal sphere. 
They find virial radii of 0.5 - 1 Mpc, the same range 
we find for for rich clusters. lYee fc EUingsoiJ l|2003f) 
examined i?2oo as a function of cluster richness Bgc (the 
amplitude of the galaxy cluster center correlation func- 
tion measured for each cluster, scaled by a luminosity 
function and spatial distribution) for 16 rich clusters with 
spectroscopically confirmed members. They measured 
-R200 by applying a singular isothermal sphere m odel to 
the v elocity dispersion data (Carlberg, Ycc & EU ingsonl 
[1993), and find i?2oo B°-^^'^^°-^^ . This scaling matches 
well with what we find, although with greater scatter 
than our result. With our techniques, we have been able 
to use photometric data to confirm this relationship for 
rich clusters and extend the relationship to much less rich 
systems. 

To compare equivalent regions of clusters of different 
richness, i?200 ™ay be used as an aperture within which 
to compare the properties of cluster members. The space 
density distribution of galaxies within i?2oo is well de- 
scribed by an NFW profile, and the derived concentration 
parameter varies with cluster richness. We examined the 
population distribution function of cluster galaxies to de- 
termine how the luminosity function of cluster members 
changes both radially and with cluster richness, using 
our determination of i?2oo to compare clusters in a wide 
range of richnesses. The radial variation of the luminos- 
ity function of cluster galaxies is similar in clusters of 
all richnesses, but does depend on the cluster richness. 
It is important to note that we can still detect a signal 
from the clusters at 2/i~^ Mpc from the cluster, even for 
poor groups. That there is still a significant over-density 
at these distances suggests caution when measuring the 
contribution of the background in an annulus centered 
on the cluster. 

We find that the central galaxies of clusters are distin- 
guishable from the rest of the cluster galaxy population, 
as has been noted by others in both observational stud- 
ies and in theoretical models. The BCG population is 
clearly evidenced in the luminosity function within R^q, 
and more dramatically in the LP of the central region 
of clusters, as a "bump" at bright magnitudes that rises 
above the LP of the other cluster members. As the rich- 
ness of the cluster increases, the BCG population be- 
comes brighter but contributes less to the overall cluster 
light, in agreement with the results of if-band observa- 
tions of BCGs in 93 X-ray selected systems studied by 
iLin fc Mohil 1)20041) . We find that the luminosity func- 
tion of low- and intermediate-richness systems cannot be 
well described by a Schechter function when the BCGs 
are included; this function provides an acceptable fit only 
for the richer systems where the fractional contribution 
of the central galaxy to the LP is small. The total cluster 
luminosity function within i?200 ^s^^ modeled by the 
sum of a Gaussian for the central galaxy and a Schechter 
function for the satellites over the whole range of Ngais- 
The trends of mean BCG luminosity and fractional con- 
tribution of the BCGs are in good agreement with the 
models of iZheng etal] l|2004|) . We find that the mean 



BCG luminosity scales with mean log Ngais as Ngais^^^, 
a similar though slightly steeper result to what is seen in 
other studies. 

The low-richness objects identified with our cluster 
finding algorithm have properties that are quite differ- 
ent from rich clusters. Systems with Ngais < 5 and Ngais 
> 7 have clearly different radial profiles and luminosity 
functions. The transition region 5 < Ngais < 7 has in- 
termediate properties. One possible contribution to this 
difference is that the lowest-Ngais objects are also the 
most likely to suffer from the effects of misidentifying the 
cluster center and/or the the brightest cluster member. 
In addition, the Ngais < 8 objects have a strong selection 
function: maxBCG is finding systems where specifically 
only a few red galaxies are within Ih^^ Mpc of the galaxy 
identified as the BCG. It is interesting that demanding 
so few galaxies of that color be in the neighborhood has 
the effect of finding objects that have very concentrated 
radial profiles, are located in isolated regions, and are 
dominated by bright galaxies. 

A possible explanation for these very low Ngais sys- 
tems is that at least some of them are fossil groups. 
Observationally, a fossil group has the X-ray lumi- 
nosity of a group or poor cluster, but in the opti- 
cal, only a highly luminous early-type galaxy without 
bright neighbor s is observed . The first such object 
was detected bv lPonman et a l. (1994); subsequent study 
has shown that there are a few other galaxies in this 
group, but all are significantly faint er (Am > 2.5mag) 
than the first-ranked galaxy ((.Tones. Ponma n fc Porbea 
i200Qj). Other fossil groups ha ve b een i dentified as 
well (El ulchaev & Zabludoff 1999; V ikhhnin et al.lll999l: 
iR.omer et al.iii2000i: .Matsushita. 200L) . Such systems are 
thought to be the end-product of groups in which most 
of the galaxies have merged, producing a highly lumi- 
nous central galaxy mostly alone in a halo of hot X-ray 
gas. This scenario is supported by the observations that 
the central galaxies emit a very high fraction of the total 
optical light of the group, and that there is a dearth 
of L* galaxies in the central region of these systems 
(Jones et al. 2003). 

At fixed halo mass, there is a theoretical expectation 
that there should be a correlation between halo formation 
time and the tota l number of galaxies a bove a given lumi- 
nosity threshold (jZentner et al J 120041) . This correlation 
may select a special population of of clusters at low Ngais ■ 
We speculate that in high mass halos, all systems have 
bright red galaxies regardless of their formation time, and 
are equally well found by the maxBCG algorithm, while 
in smaller, group-mass halos the requirement to select 
red galaxies may select only the earliest forming halos. 
These low Ngais systems would be expected to be early- 
forming halos with higher than average masses for their 
Ngais, whose outer satellites have already merged with 
the central object. This theory is consistent with the 
LP shape and the extremely steep radial profiles seen 
here, as well as with the higher masses that are indi- 
cated by cluster-mas s correlation function examined by 
ISheldonet all ll2005^ . 

Many of these systems may also be classified as com- 
pact group s (CGs), such as those studied bv iHicksonI 
[1982, 1993 and others. Compact groups have been iden- 
tified in SDSS data by Lee et aL (20041 . They investi- 
gate the local environment of CGs, and find that while 
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on average the number density of surrounding galaxies is 
comparable to the local environment around field galax- 
ies, there is considerable scatter to both more and less 
dense environments (see their Figure 9). By demand- 
ing few other bright, red galaxies nearby, the maxBCG 
algorithm may be preferentially selecting fossil groups, 
which would also look like those CGs in environments 
with a low number density of neighbors. Further study 
is needed to understand this interesting set of systems. 
More investigation can be done using SDSS data, but it 
also would be interesting to look for an X-ray signal from 
a stacked set of these low-richness objects. 

We draw attention to the range of luminosity functions 
for galaxies within R^q of clusters of different richness, 
and also to the differences in the LF seen radially. Such 
variation makes it difficult to compare between differ- 
ent cluster catalogs, where different proxies for mass and 
aperture are in use. The LF of clusters in a mass range 
are different when different fractions of i?200 sampled, 
and different again when a fixed aperture is used for each 
cluster. Nonetheless, we do see similar results to those of 
other authors, including an upturn at the faint end for 
rich clusters. This effect comes primarily from galaxies 
located in the outer regions of the clusters. 

We are currently working to compare our results with 
simulations to show how our space density-based i?2oo 
lates to the mass density-based R200 commonly used in 
N-body models. Also, we will examine how our richness 
parameter, Ngais, relates to mass M200 from the simula- 
tions. Further work will be done to use the rich SDSS 
dataset to explore the distribution of many properties of 
cluster galaxies in addition to those considered here. 

Using a small subset of the SDSS data is sufhcient to 
determine the scaling of i?200 with cluster richness, and 
to robustly measure the luminosity functions and radial 
profiles in small Ngais bins. Using the full SDSS dataset 
will allow us to study ^ 500 of the richest clusters and 
over 25,000 groups, and will allow cluster detection to 
higher redshifts. We will check for evolution in our scal- 
ing relationship by examining the i?200"Ngais trend in dif- 
ferent redshift samples. Preliminary results indicate that 
the scaling is the same, but with higher normalization for 
lower redshift objects. To compare the properties of any 
stacked set of clusters, the relationship of i?200 to Ngais is 
essential. We can now compare different mass estimates, 
including those measured from total luminosity, velocity 
dispersion, and from lensing information. We will also 
use lensing studies to compare the profiles of the lumi- 
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nous and dark components of the clusters, more closely 
investigate the possible differences in concentration of 
the two profiles, and examine the bias between luminous 
and dark matter in these dense environments. Extend- 
ing this work to include color and morphology indicators, 
we will probe many characteristics of galaxies in clusters, 
exploring the history of galaxy formation as a function 
of environment, and studying the BCG population. 

For large surveys it is typically not practical to deter- 
mine cluster mass or size using spectroscopy. Without 
assuming a model for the distribution of galaxies, we 
provide a way to determine cluster size from photomet- 
ric data alone for z < 0.3. With knowledge of how the 
size-mass scaling evolves and a better understanding of 
the scaling of -R200 with total mass, our method pro- 
vides a feasible way to measure the characteristic radii 
and masses of clusters that will be found in future large, 
high-redshift surveys. 
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